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Abstract 

We argue that regardless of the outcome of future Long Baseline experiments, 
additional information will be needed to unambiguously decide among the different 
scenarios of neutrino mixing. 

We use, for this purpose, a simple test of underground data: an asymmetry 
between downward and upward going events. Such an asymmetry, in which mat- 
ter effects can be crucial, tests electron and muon neutrino data separately and 
can be compared with the theoretical prediction without relying on any simulation 
program. 
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If the experimental evidence for neutrino oscillations and for nonvanishing neutrino 
masses consolidates further, the primary goals of leptonic electroweak physics will be 
threefold: definite identification of the flavour channels into which given initial states 
oscillate, determination of the actual number and nature of neutral leptonic states which 
participate in the mixing, and measurement of the (squared) mass differences of the 
states involved in the mixing. Obviously, the analysis of any single experiment, or subset 
of experiments, in terms of only two families can yield no more than a parametrization of 
that experiment, or group of data, and no special physical significance can be attached to 
the numerical values of the parameters. In particular, the differences of squared masses 
extracted in this way may lead to erroneous conclusions when applied to the analysis of 
other experiments. Although there is growing evidence, by now more or less accepted 
by the community, that a scenario of three flavours (e, /i, r) all of which mix strongly, is 
compatible with all data showing neutrino anomalies, a mixing pattern involving a fourth, 
sterile, neutrino Vs cannot yet be excluded. 

In the endeavour to answer these questions much hope is placed in the planned or 
forthcoming long baseline (LB) oscillation experiments M. Although these experiments 
will be very important, because of their much increased sensitivity as compared to short 
baseline oscillations, they are not sufficient to fix the parameters of the mixed neutrino 
sector. As we show in this note, LB experiments need complementary information for an 
unambiguous determination of oscillation channels and mass differences and for discrim- 
inating between the setup with three flavours and an alternative which includes a sterile 
neutrino state. We point out that the up-down asymmetry (zenith-antipode asymmetry) 
of atmospheric neutrino beams might prove to be the complementary information that 
is needed, provided the matter effects on the neutrinos coming from the antipode are 
properly taken into account. We calculate these effects by numerical integration of the 
neutrino's evolution equation and, depending on the assumed scenario and on the neu- 
trino energy, find them to be important. As the mass distribution in the earth's interior is 
well known, and as asymmetries are independent, to a large extent, of systematic errors, 
a clean analysis should be possible. 

The interpretation of the modulation of neutrino events with zenith angle reported by 
Super Kamiokande in terms of — u T oscillations, obviously, would receive support by 
observing appearance of v T in LB experiments. This could be established either directly, 
through observation of r-leptons, or indirectly, via an enhancement of the neutral to 
charged current ratio, provided no large v e appearance is found. The dependence on zenith 
angle of the atmospheric neutrinos is confirmed by recent data from Super Kamiokande. 
The analysis is performed separately on sub-GeV (p < 1.3 GeV) and multi-GeV (p > 1.3 
GeV) data samples. Electron-like [y e scattering) candidates and muon-like (z/ M scattering) 
candidates are presented separately as a function of the zenith angle (see Fig. 4 of [0). 
Looking at the electron-flavour data, one notes that there appears to be a small excess in 
the lowest cos# 2 bin. This fact allows for the Super Kamiokande data to accommodate 
some Up — > v e at low Am 2 which should be confirmed or disproved with further and more 
precise data by Super Kamiokande itself, or by the forthcoming LB neutrino experiments. 

With three generations of neutrinos, one expects a more complicated oscillation phe- 
nomenology that includes transitions between all three pairs of flavour channels. In partic- 
ular, an experiment such as Super Kamiokande which measures disappearance, might 
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in fact be seeing a combination of both — > v T and — > z/ e . In turn, if one analyzes 
the data in terms of a two-flavour scenario, say — > v T to name the most popular, one 
will extract a mixing angle and a unique, effective AM. 2 which is some convolution of the 
two physical differences AM 2 = m 2 — m| and Am 2 = m 2 — m\ but by itself has no real 
physical meaning. 

Regarding LB experiments, and keeping in mind the lesson we have learned from 
LSND 0, we already know that even if there is a large v ^ — v T appearance signal it is 
not granted that the energy distribution of the appearance signal alone would allow to 
extract Am 2 or, at least, to reduce significantly the available parameter space. To witness, 
one should keep in mind that even though LSND reports an excess of 50 events (plus a 
small background) they are not able to distinguish the low and high Am 2 cases from their 
Ufj, — u e data. 

It follows that a positive signal at any of the LB experiments cannot be taken as 
support for atmospheric oscillating into v T unless the corresponding AM 2 can be 
inferred from the same signal. The reason for this is clear: The same signal (for both 
appearance and disappearance) can be hiding a maximal — v$ oscillation with a small 
v n — v T contamination (with sin 2 (26 l ) « 10~ 3 ) or, alternatively, a three neutrino mixing 
scheme with a low — v e mixing rate. It is also important to notice that sterile neutrinos 
can only be invoked in cases where a deficit is observed, as opposed to a signal. So, sterile 
neutrinos may provide the explanation for the deficits in atmospheric or solar neutrinos, 
but LSND must indeed be — v e . At this point it may worth stressing that although 
the LSND signal is somehow under suspicion because it lies very close to the region 
already excluded by KARMEN one has to keep in mind that in their running to date, 
KARMEN sees no event indicative for — v e from the expected 3 background + 1 signal 
events (taking LSND at face value). From this result the authors conclude that not seeing 
any event allows them to exclude the LSND signal at 90 % CL. . However, had they 
observed the expected background events, this sensitivity would not have been sufficient 
to exclude the LSND result to this confidence level. 

Thus, independently of the outcome of LB neutrino experiments we are urged to look 
for complementary information if an unambiguous interpretation of their future results 
is to be possible. The complementary experiment we need must be one that allows us 
to distinguish in a clear and unambiguous manner three possibilities in interpreting the 
anomaly in atmospheric neutrino fluxes: (i) the two-flavour interpretation in terms of the 

and v T channels, (ii) the hypothesis of the channel oscillating into an otherwise sterile 
neutrino Ug, and (iii) the scenario of three strongly mixing flavours with two differences 
of squared masses AM 2 = m 2 — m 2 and Am 2 = m 2 — m\ . In the absence of such a cross 
check, none of these hypotheses can be excluded or taken for granted. 

We argue, following a proposal by Flanagan, Learned and Pakvasa [f|, that a quan- 
tity as simple as a directional asymmetry of atmospheric neutrino fluxes might prove 
to be just what we need for the purpose of distinguishing the alternatives described 
above. We feel encouraged to explore this possibility by the recent sizeable upgrading 
of Super Kamiokande which already yields better statistics than the one achieved in the 
previous, entire Kamiokande project. 
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The up-down asymmetry in e- or /i- flavour events, or Aw) ( respectively, is defined 

as 

n(/) _ f/(/) 

AV) = — — f = ue (1) 

D (f) + U (f) ' J W 

where D is the number of downward going events, produced in the atmosphere in the 
zenith, for either electron or muon neutrinos, while U is the number of upward going 
events, stemming from the antipode and having passed through the center of the earth. 
For the case of /^-flavour, D and U are 



pD(U) , 1 pD(U) 



(2) 



and similarly for electron flavour. Here N^r e \ is the initial flux of muon (electron) neutrinos 
and r = N®/N® is the expected ratio of fluxes without oscillations. This ratio varies 
somewhat with energy and zenith angle || but typical values are ~ 1/3 for Multi-GeV 
and ~ 1/1.6 for sub-GeV fluxes. 

The important point to note is that matter effects on neutrino oscillations due to the 
charged current interaction with electrons in the earth can play a considerable role. As 
was shown long ago by Mikheyev and Smirnov 0, under certain conditions the presence 
of matter can lead to a resonant amplification of the neutrino transitions, even if these 
transitions are strongly suppressed in vacuum. 

Matter effects on the oscillation of neutrino beams which cross the earth were studied 
previously in the two generation case in ||, and in the three generation case in ||. Here 
we study these effects for oscillations involving the three flavour neutrinos, u e , z/ M and u T , 
by solving the neutrino evolution equation. As the matter distribution of the earth neither 
is homogeneous nor can be modeled by means of a simple analytic function, we cannot 
rely on any simplifying assumptions and, therefore, we obtain solutions of the evolution 
equation by numerical integration. Note that with minor modifications our results can be 
extended to oscillations of anti-neutrinos, V^, and T\. 

In the presence of matter the neutrino wave function = {| u e },\ v^), \ v T )} T obeys 
the evolution equation 



with 

M{1) = U 



ij{* {v) = 7^M(l)¥^ (3) 
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m 2 2 I f/ T + I I (4) 
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Here rrii are the neutrino mass eigenvalues, and U a i is the 3x3 mixing matrix in vacuum. 
The mixing matrix relates the weak interaction states a and the mass eigenstates % in the 
leptonic sector, viz. 

3 

I v a ) = J2U ai \ ni(p,mi)) (5) 

3=1 
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with | n,i(p, 771,)) denoting the state vector of the mass eigenstate with momentum p and 
mass rrii in vacuum. It will prove to be convenient to use the following parametrization 
of the mixing matrix U in vacuum, 



U 



( 

S1C2 
\ sis 2 - 



C1C3 

- cis 2 s 3 e l 
cic 2 s 3 e lS 



c\c 2 - 
-cis 2 



S1C3 

- sis 2 s 3 e lS 
sic 2 s 3 e l 




(6) 



The background density of electrons, N e (l) induces a mass-like interaction term for 
the electron neutrino 



a(l) = 2V2G F E N e (l) 



(7) 



with E the neutrino energy. 

For a concrete application to the case at stake, these equations must be solved numer- 
ically. Note, however, that one can write down a formal solution making use of the fact 
that the mass matrix in matter, M(l), can be diagonalized at each position I by means of 
a mixing matrix Um(1)- The formal solution reads 



A (rii (I = 0) -> rij (I = L)) = P exp 



/ \ L dl \—M{l) 







2E 



where M is a diagonal matrix containing the eigenvalues of M(l) while P denotes path 
ordering. Here A (rii (1 = 0)—* rij (I = L)) is the transition amplitude for the zth mass 
eigenstate created at / = to turn into the jth mass eigenstate detected at I = L. The 
probability of one flavour eigenstate to propagate to another is then given by 



V (v a (I = 0) -> v p (I = L)) = \{v p I Um(L)A (rii (I = 0) 
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(l = L))Ul(0) 



v*)\ 2 (9) 



In the case of oscillations involving anti-neutrinos the amplitude A of finding the anti- 
neutrino (5 at I = L is obtained from an evolution equation similar to Eq.([|). The 
corresponding probabilities are obtained from Eq.@ by making the formal change U — > 
U* and a (I) — > —a(l). Note, however, that when there is a MSW resonance in the particle 
sector due to level crossing, there is no resonance in the antiparticle sector. 

In our analysis of the effects of earth on the neutrino beam from the antipode, we 
have used the density distribution as given by what is called the preliminary reference 
earth model |fT0| , Pe( t )i r denoting the distance from the center of the earth. According 
to this model, Pe(t) increases from an initial value of 1.02 gr/cm 3 in the earth surface to 
its maximum of 13.1 gr/cm 3 in the center of the earth. There are basically nine regions 
in which pe(r) varies continuously. The discontinuities of pe(r) at the borders of these 
regions are described by step functions. The most pronounced of these jumps (all of which 
are rather small) takes place at the border between the mantle and the core where Pe(t) 
changes by 4.3 gr/cm 3 . This happens at a distance r= 3480 km from the center of the 
earth. 

In our calculation we assume that the ratio of the electron density to the nucleon 
density is everywhere the same and is equal to 1/2, Ne = is the nucleon 

density and is given by = PeNa with Na Avogadro's number. 
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The result of this tedious calculation can be summarized as follows: For upgoing neutri- 
nos originating from the antipode, matter effects affect substantially neutrino transitions 
if 

10 3 GeV/eV 2 < , < 10 5 GeV/eV 2 (10) 



or equivalently 



Am 2 ' AM 2 



10-" eV < , < 10- eV . (11) 



(As Am 2 <C AM 2 we have taken mj — m 2 ~ AM 2 ). The magnitude of the interval (10) 
or (11) is determined by the range of values of the electron number density that we find 
in earth when crossing it along a diameter. Remember that matter effects are important 
only if the interaction energy squared (7) 

is comparable to or larger than either Am 2 or AM 2 , cf. Eqs. (13) and (14) below. The 
enhancement typically shows up in the dependence of the probability of a given transition 
on the neutrino energy as an irregular sequence of two or three well pronounced local 
maxima with different heights. 

For E/Am 2 (M 2 ) > 10 5 GeV/eV 2 or, equivalently, Am 2 (M 2 )/E < 10- 14 eV, the reso- 
nant densities \m 2 — m 2 \ cos(2#)/ (2\/2GpE) are much smaller than the electron density in 
the earth which, as mentioned before, varies from 1 to 13 gr/cm 3 . In this case the charged- 
current interaction with electrons in the earth dominates and oscillations are suppressed. 
For E/Am 2 (M 2 ) < 10 3 GeV/eV 2 or, equivalently, Am 2 (M 2 )/E > 10~ 12 eV, the resonant 
densities are much larger than the electron density and neutrinos oscillate like in vacuum. 

When a beam of electron neutrinos crosses the earth there can potentially be two 
resonances, in the v e — > and the v e — > v T channels. In the case of muon neutrinos there 
can only be one resonance, the one in the — > v e channel. In the parametrization (6) 
the resonance densities are given by 

N (res) Am 2 cos 29i (res) 

Ve-tfp 2\[2EGp Vn-*"e> \ ^1 

Ve Ut 2\[2EGp y ' 

These resonances would show up wonderfully if the vacuum mixing angles, sin and 
sin #3 1 were small, fLl] , provided they are sufficiently separated 0], i.e. provided 



AM 2 - Am 2 > (AM 2 sin 6 3 + Am 2 sin : 

It can be shown that the resonances take place at densities at which the differences of 
local mass eigenvalues M2(Z) — M\(l) and M^{1) — Mi(Z) are minimal. One advantage of 
the parametrization (6) we have chosen for the mixing matrix is that Mi(Z), M2(Z) and 
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M 3 (/) do not depend on the CP violating phase 5 and therefore the very existence of these 
resonances as well as the electron number densities at which they occur are independent 
of S. Indeed, the eigenvalues of M(l) are given by 

Mf (Z) = m\ + - - ^VA 2 - 3B (cose; + v^sinc;) (15) 

Mf (Z) = mj + — - ^VA 2 - 3B {cosu - v^sincj) (16) 

M 3 2 (Z) = m\ + - + \jA 2 - 35COSW (17) 

where 

A = 2Am 2 + AM 2 + a(l) (18) 

B = (AM 2 + Am 2 ) Am 2 + a(l) [(AM 2 + Am 2 ) cos 2 6 3 + 

Am 2 (cos 2 9 3 cos 2 9 1 + sin 2 3 )] (19) 

C = o(0 (AM 2 + Am 2 ) Am 2 cos 2 # 3 cos 2 #i (20) 

1 2A 3 - 9AB + 27C . , 

u; = - arccos 5 (21) 

3 2{A 2 -3BY . 

Note, however, that the transition probabilities do depend in a nontrivial way on S. 

We now turn to the analysis of the up-down asymmetry for both electron and muon 
neutrinos. We calculate this asymmetry for the three scenarios listed in the introduction, 
i.e. the two-flavour interpretation, oscillation into a sterile v$ channel, and strong mixing 
of three flavours with two rather different mass differences, all of which would give similar 
and hardly distinguishable signals in LB experiments. 

The two-flavour mixing is the simplest possible case, although it is bound to explain 
the atmospheric neutrino anomaly only, while discretely keeping silent about the solar 
neutrino deficit and the LSND result. Here, the electron neutrino flux is unaffected and 
therefore, the corresponding electron asymmetry A^ is essentially zero. In order to 
calculate the muon asymmetry, we have to take the limit sin Q\ — > and sin # 3 — > in the 
three neutrino case outlined before. Although A^ is nonzero, matter effects are negligible 
in this case. 

In the scenario of three strongly mixing flavours the larger of the two differences 
AM 2 = m 2 — m\ = .3 eV 2 is found to account for both, the atmospheric anomaly at 
low energy as well as the observations of LSND. The second mass difference Am 2 = 
m\ — m\ ~ 10~ 3 — 10~ 4 eV 2 explains the solar neutrino deficit but it is also called to play 
an important role in the atmospheric neutrino anomaly for upward going events. In this 
case, as of the two would-be resonances the one corresponding to v e — > is reached in 
the core of the earth, matter effects are dominant. 

Finally, another way to accommodate all observations is to introduce at least one 
additional neutrino state v$ which mixes with the active neutrinos but otherwise is sterile. 
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In Ref [n] the four possible types of neutrino spectra were considered. It was shown there 
that only two out of all possibilities are compatible with the existing bounds. In these 
two schemes, the four neutrino masses are grouped in two pairs of closeby masses, the 
two groups being separated by a gap of the order of 1 eV. Thus, it is assumed that 
Am^j = 1 eV 2 (we again use the abbreviation Am| = mf — m 2 ). In one scheme, say 
(A), the choices are Am^ ~ 10~ 3 eV 2 (in view of explaining the atmospheric neutrino 
anomaly) and Am| 3 ~ 10~ 10 or 10~ 5 eV 2 (needed for the suppression of the solar electron 
neutrinos). Here the numbers refer to the vacuum oscillation solution and the MSW 
solution, respectively. In the other scheme, say (B), the roles of Am^ and Am 2 3 are 
reversed. A detailed analysis of these models as well as the bounds following from them 
can be found in |I5|. 



What is important to notice is that the effective Hamiltonian of the interaction of 
neutrinos with matter in the case of models containing active and sterile neutrinos, has 
an additional neutral-current term, apart from the usual charged-current term, which is 
proportional to the background density of neutrons, and therefore has half the charged- 
current term strength. 

Cutting this long story short, in both models (A) and (B), oscillations between the 
electron neutrino and the sterile neutrino are not affected by matter effects. In the case 
of model (A) the v e — * oscillations are in the appropriate range for resonant effects to 
occur while crossing the earth and, therefore, electron neutrino beams will oscillate and 
get depleted substantially. The resulting electron asymmetry would be strikingly 
different from the three-flavour case: its property of having the opposite sign in the two 
scenarios would make it a unique and easy criterion to distinguish these alternatives. 
Furthermore, in these models with a sterile neutrino A^ e > and A^' will differ only slightly 
from one another; they are basically similar and always have the same sign. 

Figure 1 shows our results for the electron neutrino asymmetry A^ e \ Fig. 2 gives the 
corresponding results for the muon neutrino asymmetry A^\ as a function of the neutrino 
energy and for the three models described above. It is clear that given sufficient statistics 
the three scenarios can be clearly distinguished by both, the energy dependence and the 
relative signs of the asymmetries A^ and A^ e \ In particular, it is noteworthy that in 
the case of the electron asymmetry, for which the two-flavour oscillation model predicts a 
vanishing result, both observation or non-observation of such an effect would be relevant. 
Finally, although we have restricted our analysis to these three models, it will be easy 
to calculate the expected asymmetries also for other scenarios or for a different choice of 
parameters. 

In conclusion, whatever the outcome of future Long Baseline neutrino experiments 
will be, complementary information will be crucial for an unambiguous interpretation of 
their results. A complemementary experiment is needed in order to distinguish in a clean 
manner between the — > u T , —>■ v$ and the three-flavour alternatives in interpreting 
the atmospheric neutrino anomaly. For this purpose, we proposed to use the up-down 
asymmetries of electron and muon neutrinos that are detected in underground detectors 
and that we calculated by numerical integration of the evolution equation in earth. These 
asymmetries have the virtue that they do not rely on any simulations and that they require 
only identification of charged particles and measurement of their energies and directions. 
Contrary to the zenith angle dependence, the asymmetries test electron and muon neutrino 
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data separately and do not require knowledge of the fluxes. Given sufficient statistics, 
the asymmetries and Aw, when added to the results of forthcoming long baseline 
experiments, will help to distinguish the different scenarios that are being considered and 
will allow for an unambiguous determination of the mixing matrix and the differences of 
squared neutrino masses. 
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Figure 1: Muon neutrino asymmetry, versus energy (GeV) for the two flavour mixing 
scheme (long-dashed line), the three strongly mixed flavours scenario (solid line) and the 
three active plus one sterile neutrino model (A) (dotted line) and (B) (short-dashed line) 
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